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It is an intriguing possibility that dark matter (DM) could have flavor quantum numbers like the 
quarks. We propose and investigate a class of UV-complete models of this kind, in which the dark 
matter is in a scalar triplet of an SU(3) flavor symmetry, and interacts with quarks via a colored 
flavor-singlet fermionic mediator. Such mediators could be discovered at the LHC if their masses 
are ~ 1 TeV. We constrain the DM-mediator couplings using relic abundance, direct detection, and 
flavor-changing neutral-current considerations. We And that, for reasonable values of its couplings, 
scalar flavored DM can contribute significantly to the real and imaginary parts of the Bs-Bs mixing 
amplitude. We further assess the potential for such models to explain the galactic center GeV 
gamma-ray excess. 


1. INTRODUCTION 

In recent years, model builders have entertained the 
idea that dark matter (DM) comes in three generations 
like the matter particles of the standard model (SM), 
and that its interactions with the SM are governed by 
an approximate flavor symmetry. In the design of such a 
model, one must decide whether the dark matter carries 
quark or lepton flavor. In this paper we focus on quark- 
flavored dark matter, which has previously been studied 
in refs. m-m- A common element of such models is 
the presence of an additional new particle, the mediator 
that carries the quantum numbers of the standard model 
quarks, to which the dark matter couples. 

One must also decide whether the dark matter is a 
fermion or a scalar (implying the opposite choice for the 
mediator). So far, previous studies have assumed the for¬ 
mer, which we refer to as FDM. The scalar case, which 
we call SFDM, has some distinctive features that deserve 
investigation; we aim to fill this gap in the present pa¬ 
per. One difference is that the colored fermionic media¬ 
tors X bave a larger production cross section at the Large 
Hadron Collider, improving the prospects for their dis¬ 
covery or tightening constraints on their masses relative 
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to scalar mediators. 

Another difference is that scalar DM (j) can couple to 
the Higgs by the renormalizable operator A|0p|iLp that 
leads to Higgs portal interactions. We will show that 
this naturally dominates over the mediator interactions 
for setting the relic abundance and indirect signals for 
light dark matter, putting SFDM on a similar footing to 
minimal scalar dark matter in these respects. However, 
for heavy DM with mass ~ 450-1000 GeV, media¬ 
tor exchange with annihilation to tt can dominate over 
Higgs portal annihilations. Moreover, the mediator ex¬ 
changes can lead to important effects for direct detection 
and flavor-changing neutral-current (FCNC) processes. 

A further motivation for our study arises from indica¬ 
tions of an excess of multi-GeV energy gamma rays from 
the galactic center (GC), whose origin is not obviously 
tied to known astrophysical sources mHis]. There has 
been considerable interest in dark matter annihilations 
into standard model particles as a possible explanation 
of the signal, including FDM models [5] . Here we update 
the status of scalar dark matter annihilations through 
the Higgs portal to fit the GG excess, taking account of 
newer data sets provided by refs. [Min]- 

In the following we dehne the models (section]^, de¬ 
rive constraints on the mediator masses/couplings (m^ 
and Kij) from the LHG (section]^, and show the impli¬ 
cations for the couplings from requiring a thermal origin 
for the abundance (section]^. Constraints from indirect 
detection, as well as the tentative evidence for the GC 
excess, are examined in section followed by a study of 
direct detection (section 1^. Additional bounds on 
and Kij from FCNC searoies are presented in section]^ 
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In section]^ we illustrate the range of possible effects in 
this model A summary of our findings is given in the 
concluding section 


2. MODELS 

The largest quark flavor symmetry group is a prod¬ 
uct of three SU(3)’s, SU(3)QxSU(3)„xSU(3)£i, where Q 
denotes quark doublets and u,d the weak singlets. If 
we took the dark matter triplet to transform under one 
of these SU(3)’s, it would be natural to invoke minimal 
flavor violation (MFV) [TH] to suppress FCNCs in our 
model. However, this transformation choice is not neces¬ 
sary; it is more general to assume that the DM transforms 
under its own SU(3)0 group [TU], which like the others 
gets spontaneously broken by the mechanism that gener¬ 
ates the Yukawa couplings. We adopt this more general 
approach here. 

This leads to three possible models, characterized by 
the quantum numbers of the mediator particle x- All of 
them have interactions of the form 

(t>*i xPL,Rqj + h.c. , (1) 

where qj stands for quark doublet Qj or singlets uj, dj, 
and Pl.r projects onto left- or right-handed states (left 
for Qj and right for uj, dj). We will denote the models 
by Q,u,d, according to the kind of quarks which ap¬ 
pear in Q, and which the mediator must resemble in 
most respects. The differences are that the mediator has 
no generation index, and it is vector-like, having a mass 
> I TeV (see sect. below) that is independent from 
electroweak symmetry oreaking. 

In addition to the interactions with quarks, scalar dark 
matter can couple to the Higgs via 

\H\^ , ( 2 ) 

where Ay is Hermitian. At scales above that where flavor 
symmetry is broken, one expects the flavor-symmetric 
form Ay = Ao Sij , but this gets flavor-breaking radiative 
corrections that we will discuss below. 




Figure 1: Decay of the heavy DM state to the lightest one: 
(a) tree-level diagram for (^2 —> ?2fZi0i; (b) two-loop diagram 
for (f )2 —>■ 4>iX- 


2.1. DM mass spectrum and couplings 

Like the coupling ([^, we expect the mass matrix for 
scalar triplet dark matter to be flavor-conserving at high 
scales, but corrected by flavor-breaking self-energies at 
one loop, and also by the contribution from (|^ due to 
electroweak symmetry breaking: 

= ml Sij + m\ (AA^)y -h u^Ay -h ... , (3) 

where v = 174 GeV is the VEV of the complex Higgs 
field. A priori, there are no restrictions on the structure 
of Ay nor do we know whether mg dominates over the 
other contributions. For simplicity of notation we will 
henceforth take Ay to denote the matrix of couplings 
in the DM/quark mass-eigenstate basis, and allow the 
spectrum of DM states to be arbitrary, with being the 
lightest. 


2.2. Decays of excited DM states 

A priori we have three dark matter particles since (Si 
is a triplet. As long as the mediators are heavier than 
the DM, the decay </)—;► yg is forbidden and the Z 2 sym¬ 
metry under which both (jj and x ^.re charged guarantees 
the stability of (j)i. However, if there are mass splittings, 
as we generically expect there to be, then only the light¬ 
est state is stable, since a heavier one (j )2 can decay via 
(/>2 q 2 qi(Si- Even if the mass splitting is too small to 
produce the quarks, they can be virtual in a two-loop 
diagram to give (j )2 —t as shown in hg. (In fact 

the photon must be off-shell since the effective operator 
d^(j) 2 dv<t>iP^'^ vanishes for on-shell photons, but we can 
have for example (j )2 —t (Sie'^e~.) 


3. LHC CONSTRAINTS ON MEDIATORS 

The colored fermionic mediators of the model may be 
produced at the LHC, giving constraints on the mass 
of the new particle. Pair production of the mediator, 
with subsequent decay x q4>i contributes to a signal 
characterized by hnal-state jets and missing transverse 
energy, denoted $t- This is also the signature of squark 
and gluino production in the supersymmetric extension 
of the SM (SUSY). A recent ATLAS search for squarks 
and gluinos in this final state was presented in ref. [5D] . 

Signals corresponding to different jet multiplicities are 
sensitive to the production of mediators that couple ei¬ 
ther to tops or to light quarks. A mediator that couples 
to light quarks has an identical signature to the light 
squark in SUSY, namely two jets and fij,. However, a 
colored fermion has a larger production cross section than 
a scalar. The signature of a t-coupled mediator is more 
similar to that of gluino production, having a final state 
with higher jet multiplicity; the decays of t and t in the 
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all-hadronic channel result in up to six jetsd 

For this analysis, we simplify the model by ignoring the 
distinction between DM flavors (valid as long as their 
masses are much less than m^), and we allow for two 
mediators Xu and Xd- These can represent either the two 
components of the SU(2)-doublet x in the Q model, or 
the SU(2)-singlet mediators of the u oi d models. The 
interaction terms can then be written as 

Aui (l>XuPL{R)Ui + \di (t>XdPL{R)di + h.c. (4) 

We used MadGraphS m to calculate signal cross sec¬ 
tions and to generate parton-level events. Implementa¬ 
tion of the model in MadGraph is achieved with Feyn- 
Rules [22] . 

The electroweak contributions to the mediator produc¬ 
tion cross section are highly subdominant to the QCD 
process. Fig. [^a) shows the leading-order (LO) cross 
sections for the two subprocesses, verifying that the elec¬ 
troweak contributions may be neglected, as one would ex¬ 
pect. The limits on Xu and Xd are thus equally applicable 
to the mediators of the u and d models, respectively. 

We first consider the Xd mediator. Ref. [20] provides an 
upper limit on the pair-production cross section for light 
squarks as a function of their mass. Under the assump¬ 
tion that the signal topology does not differ substantially 
for the mediator signal, we calculate the cross section 
for XdXd production and translate this limit to a 95 % 
c.l. bound on the mediator mass, applying a RT-factor to 
account for higher-order corrections. The hadronic pro¬ 
duction mechanism of x is the same as for any colored 
fermion: we therefore estimate the iF-factor to be the 
same as for tt production, and obtain a value K = 1.5, 
comparing the NNLO value of the top pair-production 
cross section at = 8 TeV |23l [24] to the value cal¬ 
culated at LO. The limit is shown in fig. ib). We find 
that the mass of a down-coupled mediator is constrained 
to be > 920 GeV, regardless of the coupling strength, as 
long as X decays within the detector. 

In the case of Xu, the event topology of the signal may 
be substantially different than that from squark decays 
(other than i) due to the possible decay channel Xu ^ 

If this channel is suppressed, the signal is identical to that 
of Xd, and the same limit > 920 GeV applies. A dif¬ 
ferent approach is necessary for the Xu ^ t channel. In 
this case we use the ATLAS upper limit on the visible 
cross section, defined as the product of (cross section) 
X (reconstruction efficiency) x (signal acceptance), in 
other words, an effective cross section for the number of 
signal events observed. To obtain a limit, we simulate 
full events with hadronization and detector simulation 
in order to determine the signal acceptance and recon¬ 
struction efficiency of the mediator signal. Details are 
provided in appendix [A| 


^ Events with leptonic and semi-leptonic top-antitop decays are 
rejected in the analyses 


The resulting limit on Xu 4>'t, obtained from the 
95% c.l. upper bound on the visible signal cross section, 
is shown in fig. The exclusion region in the (m^u ’ ''^ 0 ) 
plane is shown for the signal regions (SRs) with the high¬ 
est sensitivity, and thereby the greatest exclusion reach. 
These correspond to SRs having four jets, with both loose 
and medium-level kinematic cuts (4jl and 4jm), the five- 
jet signal region (5j) in the case of 100% decays to top 
quarks, and the two-jet signal region, with medium-level 
cuts (2jm), for the other cases. For light DM, < 400 
GeV, the limit corresponds to a lower bound on a top- 
coupled mediator mass of ~ 1 TeV in the case of decays 
exclusively to tops. The bound relaxes with branching 
fraction to ^ 900 GeV; the exclusion by the 4j SR is re¬ 
laxed, while that of the 2j SR becomes stricter, as the 
branching fraction to tcj) is decreased.^ 


4. RELIC ABUNDANCE 

Since our DM candidate is a complex scalar, its particle 
and antiparticle are distinct and it could therefore be an 
example of asymmetric dark matter, whose abundance 
arises through the generation of a particle-antiparticle 
asymmetry in the early universe. However, this would 
require a more complicated model, so we will assume that 
such an asymmetry is negligible and that the relic abun¬ 
dance comes from thermal freezeout of the annihilation 
processes. These can proceed either through t-channel 
exchange of the mediator x or the (j)*(j)j\H\'^ coupling, 
as shown in fig. It will turn out that the former is 
the dominant process only in models with annihilation 
to top quarks, and with exceeding some minimum 
value to be determined. We consider x-mediated anni¬ 
hilations first, and subsequently treat the Higgs portal 
scenario, constraining An as a function of in order 
to get the observed abundance of dark matter. 


4.1. Mediator dominance 

We begin by evaluating the amplitude in fig. |^a) for 
4>i4’k —t qidj- In general, the final-state quarks could be 
different from each other, and likewise the initial dark 
matter flavors could be distinct. To simplify the kine¬ 
matics we will evaluate the cross section in the approx¬ 
imation that the DM mass splittings and quark masses 
are small compared to the average m^. At kinematic 
threshold, where the DM is at rest, the spin-summed. 


^ As the third-generation coupling is taken to zero, the limiting 
value of the lower bound on is slightly lower than in the 
light quark case, fig. We adopt the latter constraint, as the 

discrepancy is most likely a result of using different simulation 
and reconstruction methods than those of ref. [20| . as well as 
some subtler differences between the analyses. 
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Figure 2: Left (a): comparison of the relative magnitudes of electroweak and QCD contributions to mediator production cross 
section at the LHC for c.m. energy ^/s = 8 TeV. Right (b): lower bound on the mass of a colored mediator coupling to light 
quarks, resulting from the upper limit on XX production cross section in final states with jets and 








Figure 3: ATLAS constraints on DM mass versus up-type mediator mass for different branching fractions of x ^ 4’^ 
opposed to decays into light quarks): BF = 1, 0.5, 0.3, and 0.1 from left to right and top to bottom. Regions below and to the 
left of the dashed curves (envelope of exclusion from signal regions with 2, 4 and 5 jets) are excluded. 
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flavor multiplicity. If they are all degenerate, then 

must be increased by a further factor of relative to 
the complex singlet case. The exact value required will 
depend upon the mass splittings of the DM matter states 
and the thermal history. In particular, if the heavier DM 
states decay before freezeout, they will not contribute to 
the final abundance, whereas if they decay afterwards, 
they will. The range of possibilities is covered by the 
three curves shown in fig. [^b). 


Figure 4: Processes contributing to thermal freezeout of scalar 
dark matter: (a) by exchange of the mediator x\ (b) through 
Higgs portal interactions. Possible decays of h in the lower 
diagram are not shown. 


squared matrix element is 




2^2 2 

A * ,2 ^<1 (^0 - 

(m2 +m2 -m2)2 


( 5 ) 


including the sum over colors. The annihilation cross 
section is then given by 


(JVre\ 


327rm^ 


\/l - {mqlm^y 


( 6 ) 


4.3. When can mediators dominate? 

In order to determine in which cases the Higgs portal 
interactions dominate over mediator exchange for setting 
the relic density, we note that the couplings A^- can nat¬ 
urally be no smaller than typical values generated by the 
loop diagrams shown in fig. One could fine tune the 
bare value of Xij to cancel the loop contribution, but 
in the absence of such tuning one would expect a min¬ 
imum magnitude for A^, which we estimate by taking 
the leading logarithmic contribution and evaluating the 
log between the DM mass scale of order 100 GeV and 
a UV scale A = 10 TeV, which we take to be the min¬ 
imum scale of validity for our model, considered as a 
low-energy effective theory. In this case, = 9, 

and by evaluating the loop we get the estimate 


To get the right relic density, we can match this to the 
value found in ref. [2S], where the required cross section 
as a function of mass is derived. More specifically, for 
complex scalar DM, the required cross section is twice 
as large as that given in [^, where self-conjugate DM 
was assumed. Moreover, if the higher-mass DM states 
are in thermal equilibrium at the time of freeze-out, we 
must multiply the fiducial value of the cross section for 
a single complex scalar by the total number of complex 
DM components. 

The result is shown in fig. [^a), assuming a mediator 
mass of = 1 TeV, and considering only the case where 
the final state is tt, since for the lighter quarks, the sup¬ 
pression by leads to nonperturbatively large values of 
Aij. Thus mediator exchange can only be the dominant 
contribution to annihilation in the Q and u models. 


4.2. Higgs portal dominance 

In the case where the Higgs portal interactions domi¬ 
nate the dark matter annihilation cross section, the re¬ 
quired value of Xij can be deduced by rescaling with re¬ 
spect to real scalar singlet dark matter, which has been 
studied in detail in many references, including [26] . Since 
the abundance scales as l/{crv) ~ l/Af^ , we must increase 
Xij by a factor of for complex (j) relative to a real 
singlet, to compensate for the doubling of the number 
of degrees of freedom. In our model, there are actually 
three complex scalar dark matter states, because of the 
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ml, y = yd 
ml, y = yu ’ 


( 7 ) 

where Aij is the 0X9 coupling and yij the Yukawa cou¬ 
pling relevant to the particular model of interest; = 174 
GeV is the complex Higgs VEV. Which Yukawa matrix 
appears depends upon the mediator. If the mediator is 
u- or d-like, then y = yu or y^ respectively. But if it is 
the doublet (Q-like), then we must sum over both pos¬ 
sibilities, in which case yu dominates. In either case, 
working in the basis of diagonal Yukawa matrices gives 
the approximation shown in 0. 

By substituting the value of |Aij| shown in fig. [^a) 
into eq. 0, and comparing to the value of Xij shown in 
fig- ib), we can determine when it would be inconsis¬ 
tent to assume that mediator exchange dominates over 
Higgs portal interactions. This comparison is shown in 
fig-[!l for models where rather than y^ appears in the 
loop (otherwise, the solid curve is lower by a factor of 
(mb/rrit)'^, giving no useful constraint). We assumed for 
these curves that only one DM flavor is in equilibrium; for 
higher numbers, both curves scale upward by the same 
factor, so that the values of where they intersect do 
not change. 

It is interesting to notice that the same model-building 
choices that would suppress the loop contribution 0 also 
suppress the mediator contribution to annihilation. In 
particular, the models for which y = yd in 0 are those 
where the mediator is d-like, but these have cross sections 
for 00 —)■ qq suppressed at least by ml in eq. 0 , making 
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Figure 5: Left (a): value of DM-mediator-quark coupling |Aij| needed for thermal relic abundance from annihilations via 
f-channel mediator exchange, as a function of DM mass m^, assuming (f>* (f) ^ tt and = 1 TeV (lower curves) or = 2 
TeV (upper curves). Different curves show the dependence upon how many flavors of DM are in equilibrium at the time of 
freezeout. Right (b): value of DM-Higgs cross coupling needed for thermal relic abundance from Higgs portal annihilations, 
again showing the dependence on number of DM flavors in equilibrium at freezeout. 


H 



H* 


q 

(a) 



(b) 


Figure 6: One-loop contribution to interaction, 

where q (Q) stands for electroweak singlet (doublet) quarks, 
and the routing of weak isospin is shown for (a) singlet and 
(b) doublet mediators, respectively. 


it impossible to satisfy the relic density constraint with 
reasonable values of . 

The upshot of this analysis is that only in the Q and 
u models with 300 GeV < and = 0.5- 

1 TeV can we consistently assume mediator dominance 
of the annihilation cross section. Here we have taken 
advantage of the fact that our LHC constraint on is 
weaker for ~ 300 GeV than for lighter see fig. 

(upper left). Moreover, the tree-level value of Ay can 
exceed the minimum coming from the loop estimate in 
eq. 0; thus Higgs portal dominance is always a logical 
possibility, even when not a necessity. 


5. INDIRECT DETECTION 

Annihilation of DM in our galaxy or neighboring ones 
can produce gamma rays from the decays of final-state 
particles. The Fermi Large Area Telescope (LAT) con¬ 
tinues to improve constraints on dark matter annihila¬ 
tion into various final states, from observations of dwarf 
spheroidal galaxies that are relatively uncontaminated 



Figure 7: Comparison of the value of |Aij | needed for cor¬ 
rect relic density via Higgs portal annihilations (dashed curve) 
with the value coming from loop contributions, eq. Q, in the 
case where |Ay | is large enough for mediator exchange to give 
the right relic density (also assuming that rather than yd 
is the Yukawa coupling matrix appearing in the loop). Mod¬ 
els where the solid curve is higher than the dashed one have 
annihilations dominated by the Higgs portal, in the absence 
of fine tuning. 


by baryonic background signals [27]. The constraints 
are strongest for light dark matter, whose relic density 
is higher. They are therefore relevant in the region of 
parameter space where annihilation is primarily through 
the Higgs portal. 

For rriff, below the WjZ threshold, annihilation is al¬ 
most exclusively into bb. We reproduce the Fermi limit 
from ref. m on the annihilation cross section into bb in 
%• ia), where it is relaxed by a factor of 2 due to the 
dark matter not being self-conjugate in our model. The 
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Figure 8: Left (a): cross section for (j)*4> to annihilate through the Higgs portal, for relic density, galactic center gamma ray 
excess (1, 2 and 3cr contours for three data sets: Daylan et al, CCW and Fermi), and Fermi/LAT upper limit from dwarf 
spheroidal galaxies assuming bb or WW final states. Right (b): branching fractions for (j>*(j> to annihilate into SM final states 
through the Higgs portal (adapted from ref. |26|'). 


value needed for the observed relic density (also increased 
by the factor of 2) is also shown, suggesting that masses 
below 100 GeV are ruled out. However, for rrifj, > 70 
GeV, the dominant annihilation channel is no longer bb 
but rather WW or WW* where one of the TV’s is off 
shell; fig. |^b) shows the branching fractions into differ¬ 
ent final states. The constraint on the WW + WW* 
channel is weaker by a factor 1.3, which can allow for 
somewhat lighter dark matter (m^ ^ 80 GeV) to be con¬ 
sistent with both relic density and indirect constraints. 
The actual constraint on Higgs portal models (not deter¬ 
mined by ref. [27]) interpolates between the WW and bb 
curves in the region = 70-80 GeV. 


5.1. Galactic center y-ray excess 

A possible signal in Fermi/LAT data for dark matter 
annihilation in the galactic center has been discussed by 
several groups, most recently in refs. [T5HT7] (referred to 
here as Daylan et ai, CCW and Fermi respectively). The 
Fermi collaboration itself presented preliminary evidence 
for gamma rays in excess of those attributable to known 
astrophysical sources in the central 15° x 15° region of 
the galaxy [ni nHi- Recently, new evidence has been 
presented in favor of unresolved millisecond pulsars as a 
likely astrophysical source |28H3T|, but pending a defi¬ 
nite resolution, it is interesting to explore whether dark 
matter models can consistently explain the observations. 

Here we have used a similar methodology as in ref. [32] 
to fit 7710 and its annihilation cross section {av) to the 
excess signal as characterized respectively by Daylan et 
ai, CCW and Fermi. To generate the predicted signal, 


we compute the photon spectrum from annihilation into 
SM states with the branching ratios shown in fig. [^b) , 
using spectra from ref. |33j . These are compared to the 
data to compute statistics for which the 1, 2 and 3cr- 
allowed regions are shown in fig. [^a) . 

The three data sets are not fully consistent with one 
another, and they conflict with the Fermi dwarf con¬ 
straint except for part of the CCW 3a region. This region 
however has too small an annihilation cross section with 
respect to that needed for the relic density, by a factor 
of = 1.2, which would lead to a 20% increase in the dark 
matter abundance. The experimental error in the ob¬ 
served abundance as determined by Planck is about 4%. 
A consistent interpretation would require that the actual 
excess signal be somewhat lower in intensity, as may be 
the case if part of it is due to millisecond pulsars. 


6. DIRECT DETECTION 

Our dark matter candidate can scatter elastically with 
quarks through mediator exchange, fig. |^a), leading to 
DM-nucleon scattering that is constrained by direct de¬ 
tection experiments. This can occur either by the cou¬ 
pling of (j) to valence quarks, or that to heavy quarks 
via the loop diagram fig. [^b) that enables photon ex¬ 
change. In addition, the Higgs portal coupling allows for 
—>■ (j>N scattering by Higgs boson exchange. For DM 
masses rriip > 6 GeV, the strongest current limits come 
from the LUX experiment [34]. We will first derive con¬ 
straints on the different kinds of interactions assuming 
that they do not interfere with each other significantly. 
In section [6.3| we will consider the possibility of destruc- 




































tive interference that could weaken direct detection con¬ 
straints sufficiently to allow for the indirect signals we 
discussed in section [H 


6.1. Mediator-induced interactions 


We first consider the nonelectromagnetic mediator- 
induced interaction, ft is straightforward to show that 
fig. [^a) leads to an effective operator^ 

, ( 8 ) 
ml 

where the sum over doublet components is taken in the 
Q model. When taking matrix elements of this operator 
between nucleon states, the only nonvanishing contribu¬ 
tions are from the valence quarks i = u,d, giving the 
cross section [3S] 


— 


J p,n 

47rm^ 


(9) 


for scattering of (j) on protons or neutrons, where fj, is 
the 0 -nucleon reduced mass and fp^n = up 

to isospin-related factors of order unity, depending upon 
which DM model we are considering.The LUX con¬ 
straint on these couplings is shown in fig. | 10 | [a) for me¬ 
diator mass = 1 TeV. The limit on A 13 is orders of 
magnitude smaller than values of interest for the relic 
density for the coupling to top quarks. There must be a 
large generational hierarchy in the couplings Aii, at least 
between the first and higher generations. 

Next we consider the contribution from the penguin 
diagram, fig.j^b). The loop leads to the effective photon- 
DM interaction 



Figure 9: SFDM contributions to </)-nucleon scattering via 
the bx? interaction: (a) tree-level mediator exchange, and 
(b) penguin diagram for DM-proton interaction. 


® Here we ignore contributions suppressed by niq that are irrele¬ 
vant for (^-nucleon scattering. 

^ For fp these factors are (1, 1/2, 3/2) for the u, d, Q models, while 
for fn they are (1/2, 1, 3/2), respectively. 


For the three models (u, d, Q), k is approximately given 

by 


~ 1 

IdTT^ 


QgEi lAnpln , q = uOTd 

, Q model 


in the limit ^ where Qq = ‘Ij'i or —1/3 is the 
charge of the quark and Wq. is its mass. For larger m^, 
the loop integral depends differently upon and the 
logarithm gets replaced by 



In 



ml 


( 12 ) 


where we define Ci = Wq./m^, D(x) = 

X + ei(l — x) — e 0 x(l — x), D'{x) = D{1 — x) and 


I = 


+ ( 1 -^)" 


/ ei + e^x"^ 

1 -1- e^x^A 

V 2iA2 

2D'2 ) 


(13) 


The large logarithm comes from 1/D as a: —> 0. 

The resulting photon-mediated DM-proton scattering 
cross section is given by 


IGtt 

{ml - ml) 


2\2 ’ 


(14) 


where ^ is the (/-proton reduced mass. The limits 
from this process are much weaker than those from the 
nonelectromagnetic coupling. Also, whereas that one 
bounded only An, this one applies to Aii for all the 
quark generations. Hence we take i = 2,3, (recall that 
1 = 3 represents the couplings relevant for the relic den¬ 
sity in section 4.1). Ignoring possible interference be¬ 


tween different generations, we obtain the limits shown 
in fig. [To|(b) , with solid (dashed) curves corresponding to 
i = 3(2b These couplings are somewhat weaker than the 
values leading to the right relic density in fig. [^a). 


6.2. Higgs portal interaction 


For the Higgs portal coupling, the effective DM- 
nucleon scattering cross section is given in ref. |26j : 


47r m'jpml ’ 


(15) 


where /tv = 0.303 [55] is related to the Higgs-nucleon 
coupling, and mH = 125 GeV. The LUX data can be 
used to put limits on the the coupling An, as shown 
in Fig. El The value of An needed for the observed 
relic density is also plotted, for the case where only c/i is 
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Figure 10: Left (a): LUX constraints on the coupling of dark matter to light quarks u,d, assuming mediator mass = 1 
TeV. Right (b): LUX constraints on couplings between A 13 (solid curves) and A 12 (dashed curves), for different choices of the 
quark appearing in the loop of fig. Hb), depending upon the choice of model (d, u or Q). The mediator mass is assumed to be 

■mx = L 



(GeV) 


Figure 11: Solid curves: LUX and relic density constraints 
on the FTiggs portal coupling An, taking dark matter to be 
asymmetric and accounting for accidental cancellation at the 
level of 100, 50, 25 and 10% by the mediator contribution to 
the scattering amplitude. 100 % means no cancellation, giving 
the usual exclusion curve. Dashed curve: value of An that 
would give thermal relic abundance. 


relevant during freeze-out. The LUX limit must lie below 
the relic density curve if the DM has a thermal origin. In 
addition to the allowed range > 150 GeV, there is a 
narrow window of lower masses for which the relic density 
is not exceeded around = mhl2, corresponding to 
resonantly-enhanced annihilations. 


6.3. Higgs-mediator interference 


Fig. 11 shows that the interesting DM mass range for 
indirect detection (fig. is ostensibly ruled out by the 
LUX limit. However, we have not yet taken into ac¬ 
count the possibility of destructive interference between 


different contributions to the (()-nucleon scattering am¬ 
plitude. This is clearly possible for either dark matter 
particles or their antiparticles, since the effective opera¬ 
tor ([^ changes sign under charge conjugation of (j) while 
the amplitude from Higgs exchange does not. Since eqs. 

and ( [T^ have the same dependence upon m^, it is 
particularly simple to combine them taking account of 
interference: 


Tod 


Cn = ( Ah/at—^ ± 


|Ai 






(16) 


If interference is destructive for (j) it will be constructive 
for (j)*. Therefore to have a net reduction, it is necessary 
to consider asymmetric dark matter where the antiparti¬ 
cle abundance is suppressed [35] ■ The suppression factor 
has been computed as a function r(A/Ao) = /n^ (the 

ratio of anti-DM to DM) in ref. [37] , where Aq denotes the 
value of the coupling that would give rise to the correct 
thermal relic abundance.® If the amplitude for scattering 
of DM on nucleons is reduced by the factor (1 — e), and 
that for anti-DM is increased by (1 -|- e), and the nominal 
bound is Alux, then the relaxed bound on the coupling 
is given by 

A., =-^ (17) 

[(l-e)7(l-r/2) + (r/2)(l+e)2|‘'“ 

where r = r(Aeff/Ao). Eq. 0 gives only an implicit 
definition of AeS, but it can be solved numerically by 
iteration. 

In fig. we show the modified upper limits on |Aii| 
that result from allowing for accidental cancellations that 


® Denoting x = (A/Aq)^, we are able to fit the numerical result of 
ref. 1371 to the function — logjg r = {Aqx + Ai)/(1 + A 2 X^^), 
where Aq = 0.8327, Ai = -0.8258, A 2 = -0.8737, A 3 = 
—0.8213, which is valid for x> 1. For x < 1, r = 1. 
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reduce the amplitude to 75%, 50% and 10% of its mag¬ 
nitude in the absence of the mediator contribution. It is 
clear that the range of allowed masses can be consider¬ 
ably widened relative to the thermal abundance scenario. 

The mediator diagram can have a significant effect only 
for ^-nucleon scattering, and not (j)(f>* annihilation, be¬ 
cause of the quark vector current in the effective interac¬ 
tion (|^. Its matrix element for (fxj)* —>■ qq is suppressed 
by rriq = mb for the mass range of interest, while that for 
(/)7V —> (pN suffers from no such kinematic reduction. 


7. DM-INDUCED FLAVOR EFFECTS 

We now turn to the implications of scalar flavored 
dark matter for particle-physics phenomenology, includ¬ 
ing FCNC processes, rare decays, and CP violation. We 
recall our choice of the underlying flavor symmetry group 
as SU(3)0xSU(3)QxSU(3)„xSU(3)d [TS], which is bro¬ 
ken by the SM Yukawa couplings and our new couplings 
Aij. Because the mediator y is forced to be heavy by 
LHC constraints, we do not need to rely upon a more re¬ 
strictive flavor structure such as MFV [TB] to keep flavor- 
changing neutral currents under control, as will become 
clear in this section. However, we will demonstrate the 
potential of the model to give rise to observable low- 
energy effects for values of the couplings Aij that are 
consistent with the constraints obtained in the previous 
sections. 


7.1. Flavor-changing meson oscillations 


We briefly review the formalism for AF = 2 flavor- 
changing oscillations of neutral mesons. To be concrete, 
we illustrate this for the case of AB — 2 meson mixing. 
In the basis, the mixing is described by the 2x2 

matrix M — |r, in which the mass (M) and decay (F) 
matrices are Hermitian. The physical states are [35] 

\Bl) = p\B^)+q\&) , 

\Bh) = p\B^) - q\&) , (18) 

with eigenvalues (L = “light,” H = “heavy”) 


Pl,h = Ml,h — , (19) 

in which u and F h denote the masses and decay 
widths of Bl h- In addition. 


^=± 


/ Mi*2-zrV2 V/^ 

\Mi2 — zFi2/2 


( 20 ) 


It is a good approximation to take |Fi 2 | ^ |A7i2|; then 


AM = Mh - Ml = 2\Mi2\. (21) 


In our model, the matrix element M 12 receives new con¬ 
tributions beyond the SM from box diagrams with p and 
X in the loop. 


7.1.1. B^-Bs mixing 


The DM-induced contributions to B^-B^ mixing from 
box diagrams with cp and x in the loop can be described 
by the effective operator m 


(AU)L 

1287r2 


(brPL,Rsf 


( 22 ) 


in the Q- (Pl) and d-type (Pr) models, where we used 
the approximation m^,. <C m^.^ The corresponding mass 
splitting is {AMs)dm = l(A'''A)gjTOB^/|y(1927r2 ). 

The measured value is (AMs)exp. = (11.69±0.02) x 10“® 
MeV, while the SM prediction is (AMs)sm = (11.4 ± 
1.7) X 10“® MeV |39|. These quantities are related via 


(AM,) 

exp. — 2 |(Mi2)sM + (A7i2)dm| 


(23) 


To obtain constraints on the DM contribution to 
one has to take into account a possible phase difference 
between (Mi 2 )sm and (Mi 2 )dm- But this phase differ¬ 
ence will also manifest itself in q/p, eq. (20). A rigor¬ 


ous analysis would require a simultaneous fit to the mea¬ 
sured values of AM, and arg(q/p), which is beyond the 
scope of this paper. Instead, to estimate the allowed 
size of the new contribution, we neglect any phase differ¬ 
ence, leading to |AMs|dm = (0-3 ± 1.7) x 10“® MeV, or 
|AM,|dm < 5.4 X 10“® MeV at 3a. For = 1 TeV, this 
corresponds to the limit |(AlA)f,,| < 0.19. This is smaller 
than the direct detection bounds on second-generation 
couplings shown in fig. [Toj(b). 

As can be seen from the above values of (AMs)exp. 
and (AMs)sM) the measurement of B^-Bs mixing is con¬ 
sistent with the SM prediction. On the other hand, the 
theoretical error on this prediction is sizeable, leaving 
ample room for a new-physics contribution to AM,. In¬ 
deed, if I AM, I dm saturates its upper limit, it will be of 
the same order as |(AM,)gM|. We therefore see that fla¬ 
vored DM could contribute significantly to 5,-5, mixing 
with reasonable values of the couplings. 


7.1.2. K°-K°, 5°-5°, 5S-5S mixing 


A similar analysis can be done for oscillations of the 
other neutral meson systems, K^-K°, 5°-5°, B^-B^. 
Constraints on the coefficients Cij of the effective oper¬ 
ator A“^(gi7^5iqi)^ (where A is the new physics scale) 
have been compiled for example in ref. [40] . These can 


be related to the prediction (22), with appropriate sub¬ 


stitution of quark flavors. The results are shown in ta¬ 
ble |lj For K^-K^ and 5°-5° mixing, we obtain separate 


This follows from eq. (4.1) of ref. m, accounting for the loop 
now being dominated by momenta of order and ignoring 

corrections of order 
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Figure 12: Diagrams for (a) t —>■ <l>4>*c, (b) f —>■ 07 or other gauge bosons, (c) t —>■ he, and (d) the function Firnrj,) determining 
the partial width for t —>■ 4>(l>*c relative to its maximum value. 


d' 

Re((AAt)F] 

Im((AAt)?,l 

ds 

1.1 X 10"® 

4.3 X 10"® 

uc 

7.1 X lO"'^ 

1.3 X lO"'^ 

bd 

3.6 X 

10 "^ 

bs 

8.3 X 

10 "® 


Table 1: Bounds on FCNC matrix elements of (AA^)ij with 
i ^ j from neutral meson mixing, assuming mediator mass 
= 1 TeV. Values for the first two rows are inferred using 
constraints reported in ref. [40|. For Bd,s (last two rows) we 
constrain only the modulus |(AA^)fj| using updated experi¬ 
mental and SM fit numbers from El- 


constraints on the real and imaginary parts of (AAf)| 

10 


)« 

{ij = ds,uc). For B^-B^ and B^-B^ mixing, the con¬ 
straints are given only for |(AA^)fj| {ij = bd,bs). The 
imaginary parts of (AAf)f^ can lead to CP-violating ef¬ 


fects, as we will discuss in section 7.4 


but because of the more stringent constraint on first- 
generation couplings from direct detection, this is ex¬ 
pected to be subdominant. With large couplings A ~ 3 
and light DM with ^ 30GeV, eq. ([24) would lead 
to a branching ratio of 3 x 10“®. Recent studies of this 
process in other models with flavor-changing scalar DM 
coupling to the top estimate that LHC searches could 
ultimately be sensitive to such a small branching ratio 
[45M5] . Although our choice of is ruled out by di¬ 
rect detection for a thermally produced WIMP, since (j) 
has a conserved particle number, there could be a DM 
asymmetry allowing for sufficiently small coupling to the 
Higgs for consistency with direct searches.^ 

A second class of decays ist ^ c + gauge boson, shown 
in fig. [T^b). ATLAS obtains upper limits of 1.7 x 10“^ 
on the branching ratio for t ^ eg and 4 x 10“® for t —>■ 
ug [Ml ETj. Writing the NP contribution to the t ^ c 
chromomagnetic dipole moment as 


gs Htcg {iC) G% + h.c. , (25) 


the limit on the branching ratio corresponds to ntcg < 
7.2. Flavor-changing top quark decays 1.3 x 10“^ TeV“^. In our model 


SFDM allows for a variety of rare FCNC decays, in¬ 
cluding t —>■ c{u)4>4> (if 7710 is sufficiently small), b —)■ 57, 
t —)■ {Z, h, g, j)c, and {h, Z) —>■ bs. A summary of the Fer- 
milab and LHC constraints on these processes is given in 
ref. |M]- With the exception of t ccjxj}, these are un- 
observably small, despite having no symmetry (MFV) to 
suppress them. This is a consequence of the chiral struc¬ 
ture of the interaction ([^, which causes all amplitudes 
to be suppressed by 1/r?^ and not just l/m^- 

If ^ Wt/2, the tree-level processes t —>■ c(j)(j) or 
t — >■ ucjjfj) are allowed (fig. [T^a)). For the ccjxj) final state, 
the partial width is 


ST ^ 


(AtA)„(AtA),,m? 


2xlO-^GeY-\A\l\A\lF{m^) 


4096 TT^ 771^ 


(24) 


Ktes - (A^A)t, , (26) 

implying the weak constraint (Af A)tc < 40. 

For the electromagnetic FCNC t ^ decays, CMS 
finds a limit of 1.6 x 10“^ (1.8 x 10“^ for t —>■ cy) [48j . 
This corresponds to a limit on the magnetic moment co¬ 
efficient Ktu-y 

Y Ktuy {i(J^''PL,Ru) Ff,„ + h.c. (27) 

of Ktu-y < 0.16TeV“^, and a correspondingly weaker 
limit of (AfA)t„ < 580. The best limit on t —)■ qZ also 
comes from CMS, with an upper bound of HR < 5 x 10“^ 
[45] . leading to (Af A)tc + (Af A)t„ < 785 for models with 
Q-like mediators, and somewhat less stringent for u-like. 


where the dependence on is shown in fig. H^d) and 

for the numerical estimate we took = 900 CeV. The 7 -pjjjg would also require some fine tuning of the loop contributions 

analogous formula with C —^ U applies for t — ^ ucjxj), to Xjj, according to the considerations of section|4.3| 
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The decay mode t ch shown in fig. ^c) has a par¬ 
tial width of order 

/ r;(AUA)te \^ mt(m^ - ml) 

\ J IdTT 

2 X 10-^(AUA)?,GeV , (28) 

which is far below the current sensitivity of JT < 1 GeV 
[50] for reasonable value of the couplings. 


‘I’l 



Figure 13: Diagram giving down-quark electric dipole mo¬ 
ment. 



7.3. Flavor-changing b decays 

The radiative flavor-changing processes 6 —>■ sy and 
b se+e- are described by the effective operators 

PRb)F^^ , 

09 = , 

Oio = j^{s7^PLb){iri5£) (29) 

and their chirality-flipped counterparts, O'^, Og and O'lg, 
obtained by taking Pr o Pr. The operators © 7 , 9,10 are 
induced by the Q model, while the d model generates 
O' 

'^7,9,10' 


7.3.1. b^s'y 


The decay 6 —sy has reduced sensitivity because of 
both loop and chiral suppression of the induced magnetic 
moment operator. In the Q model, to leading o rder in 
1 /m^, it is given by the diagram analogous to fig. 12 ’b). 


Qh(AU)6, ^ ^VtbVZ Or (30) 


12 


V2 


where Qb = —1/3 is the charge of the b quark, and C™ 
denotes the new DM contribution. 

In ref. |5T] , a global analysis of B decay processes was 
performed, motivated by tensions with the SM predic¬ 
tions revealed by LHGb measurements (55] . There it was 
shown that a nonzero contribution from new physics is 
preferred at Icr for © 7 , namely © 7 “ € [—0.05, —0.01]. For 
= 1 TeV, 0.5 < (AU),,* < 2.4 is the corresponding 
allowed range of couplings. 

A similar analysis can be performed for the Wilson 
coefficient ©7 which appears in the d model, where the 
1(7 range is © 7 °“ S [—0.04, 0.02] This corresponds to the 
range of couplings —0.9 < (A^A);,,, < 1.9. From both 
©^'^ and © 7 ^^^ the constraints on (AlA)bs are much 
weaker than the limit |(AlA)bs| < 0.19 we obtained above 
from Bs-Bg mixing. Hence if one imposes the Bg-Bg 
mixing constraint, our DM model cannot generate large 
enough contributions to H —>• Af^y decays to address the 
current (weak) hints of deviations from the SM predic¬ 
tions. 


7.3.2. b^s(.+r 


The new DM interactions contribute to 6 —>■ 
through b —>• sy*(—>■ or b sZ*{-^ £~^£~) (i.e., 

the y or Z is off-shell). All three of the operators in (29) 
(or their chirally-flipped counterparts) can be relevant. 

There has been a great deal of activity, both theoret¬ 
ical and experimental, concerning B fj,~ de¬ 

cays; see ref. m for a recent review. At present, there 
is a hint of new physics in ©g: at Icr, it is found that 
©^^ G [— 1 . 6 , —0.9], and remains nonzero even at 3 ct [5T] . 


Within the Q (d) model, we find that the b 
tribution to ©9 (©g) is 


sy con- 


©DM 


7V2Qb{A^A)bg 

lUGFml\Vtb\\Vtg\ ■ 


(31) 


at leading order in 1/w^. For = 1 TeV, the range 
of couplings (AlA);,s G [18, 33] corresponds to the Icr 
range of ©g*”. Similarly to the case of 6 —)■ sy, this is two 
orders of magnitude larger than the constraint from Bg- 
Bg mixing; hence the DM contribution cannot explain 
the discrepancy in ©g. For Cg™ the Itr allowed range is 
[— 0 . 2 , 0 . 8 ], again corresponding to constraints on (AlA)^^ 
that are much weaker than those from Bg-Bg mixing. 

Similar conclusions hold for all b st^t~ and b sqq 
operators. The DM contribution to the Wilson coeffi¬ 
cients is suppressed relative to the SM by 0{M'^/m^ ~ 
1%. It cannot be compensated by large values of (AlA)bs, 
due to the constraint from Bg-Bg mixing. 


7.4. CP violation 

The couplings A^ in our model can be complex, lead¬ 
ing to new sources of GP violation. They can have ob¬ 
servable effects through meson mixing in B^ decays, and 
possibly also through the electric dipole moment of the 
neutron. 


7 . 4 . 1 . Indirect CP violation in decays 


In sec. 17.11 we showed that the DM-induced contribu¬ 
tion to Bg-Bg mixing may be signihcant for reas onable 
values of the couplings (AlA)j^ ^ 0.1 in eq. (22). The 
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imaginary part of (A^ is a new source of CP violation, 
entering in g/p, eq. ( p 0 | ), 

g ^ ^ (-^ 12 )sM + (^ 12 )dm \ (- 32 ) 

P V (Afi2)sM + (Afi2)DM/ 

which is approximately a pure phase, \q/p\ = 1. In the 
SM, q/p = V:,VtslVt,V*,. 

The phase arg{q/p) can be measured using indirect CP 
violation in Bs decays. The main experimental focus has 
been on , which is the phase as measured via indirect 
CP asymmetries in Bs decays with b —>■ ccs {Bs —t 
J/'tljK~^K~, J/V'TT+TT”, DfDj). Its predicted value is 
= aTg{q/p) = — 0 . 0363 tQ;gQj 4 in the SM, while the 
measured value is —(0.015 ± 0.035) [54j. Although these 
are consistent with one another, there is ample room for a 
new-physics contribution. Given that (Afi 2 )DM 
the same order as (Afi 2 )gM! our model could contribute 
significantly to 

7.4-2. CP-violation in K^-K^, D°-D°, B^-B^ mixing 

In contrast to the B^ system, the CP phase relevant for 
K°-K° and D°-D° mixing is sufficiently well-measured 
to provide a separate constraint on the imaginary part of 
(AlA)ij for the off-diagonal elements. The upper limits 
are given in table|^ For B^-B'^ mixing, we do not present 
a constraint on the imaginary part of (AlA).y . However, 
its modulus is reasonably well-constrained, so that its 
imaginary piece cannot be too large. We therefore do 
not expect significant DM-induced contributions to CP 
violation in K^-K^, D^-D^, or B^-B^ mixing. 

7 . 4 . 3 . Electric dipole moments 

The new phases can also produce quark electric dipole 
moments through two-loop graphs like that shown in fig. 
[T^ The extra loop with W exchange is needed to get the 
products (AAl)i 4 with j ^ 1, since there are no phases 
in (AAl)ii. Because of chiral suppression, the resulting 
quark EDM is small, 

^ Im[(AAl')i 2 ] eg| Vcd mg 
^ "" {l&TT^m^Y 

= 3 X 10“^® Im[(AA'^)i 2 ] e • cm. (33) 

Given the stringent constraints on (AAl)uc and (AAl)ds 
from D^-D^ and K^-K^ mixing, this is negligible com¬ 
pared to the current sensitivity through the neutron 
EDM, 3 X 10-2® e-cm. 

8 . BENCHMARK MODELS 

Rather than trying to combine the constraints we have 
discussed to obtain allowed regions in parameter space. 



(GeV) 

All 

Aij 

comment 

1 

60 

0.01 

- 

asymmetric dark matter 

2 

63 

0.016 

- 

GC excess 

3 

100 

-0.12 

All = 0.04 

direct detection interference 

4 

200 

0.08 

- 

thermal relic/Higgs portal 

5 

700 

0.27 

00 

d 

II 

CO 

thermal relic/mediator 


Table II: Parameter values for benchmark models, assuming 
= 1 TeV. Dashes indicate where Kij can take a range of 
values. 


since we have many parameters, here we will instead give 
a few examples of allowed parameter values that illustrate 
the different qualitative possibilities of the model. For 
simplicity, we fix the mediator mass = 1 TeV, close 
to the lower limit from LHG searches. The benchmark 
models are summarized in table El 

Model 1 underscores the fact that if the annihilation 
cross section exceeds that needed for the thermal relic 
density, then we could appeal to the complex nature of 
SFDM to assume that it has an asymmetry accounting 
for its abundance. The example chosen here has An close 
to the upper limit from direct detection searches, making 
such a model close to discovery. 

Model 2 is the largest below 200 GeV allowed by 
direct detection and the thermal relic value of An shown 
in fig. [m It is marginally consistent with the galactic 
center gamma ray excess, fig. though in conflict with 
the Fermi dwarf spheroidal constraints. 

Model 3 illustrates a case that would be ruled out for 
thermal DM but is allowed for asymmetric DM due to 
destructive interference between Higgs and mediator ex¬ 
change contributions to DM-nucleon scattering. 

Model 4 is an example where the thermal relic density 
arises from Higgs portal interactions, and An is close 
to the LUX limit in fig. El again making this model 
detectable by the next improvement in sensitivity. 

Model 5 is chosen to illustrate the window of couplings 
shown in fig. where annihilation of dark matter by me¬ 
diator exchange can dominate over the Higgs portal cou¬ 
pling, without fine tuning of parameters. At this mass, 
the LUX limit upper An < 0.5 is satisfied. Moreover 
|Ai 3 I is below the direct detection limit ~ 2 shown in fig. 
[l0[b). This model could be discovered with a factor of 
3 improvement in direct detection sensitivity, via Higgs 
exchange interaction. 

The sensitivity of tests from flavor-changing particle 
physics processes generally depends upon different pa¬ 
rameters than the astrophysical ones considered above. 
For example, an observable contribution to the B^-B^ 
mixing amplitude would arise from a choice of couplings 
such that |A 2 sA 2 t,| ~ 1 in the d model. These couplings 
are unconstrained by the previous considerations. On the 
other hand it is also possible to saturate the B^-B^ mix¬ 
ing bounds using |Ai£;Ai;,| ^ 1 in the same model, but 
the constraints from direction detection shown in fig. 
rule out this possibility. This illustrates that there can 
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be some interplay between the particle physics and astro- 
physical constraints, but that in general there is freedom 
to separate them. 

The nonvanishing values of A.^ required in models 3 
and 5 do not have direct implications for meson mixing; 
rather they imply, through table |lj constraints on neigh¬ 
boring matrix elements. For example |Ai 3 | = 0.8 in the 
Q model implies |Aii| < 0.02 to satisfy mixing 

constraints. Even if An = 0 at tree level, a one-loop 
vertex correction of order 

Mil ~ -- 10”^ (34) 

is induced by exchange. Therefore no fine-tuning is 
needed to satisfy this constraint. Similarly |Ai 2 | = 0.04 
in the Q model requires |Aii| < 0.05 to satisfy K^-K^ 
mixing constraints; this is also easily compatible with the 
maximum size of loop corrections. 

9. CONCLUSIONS 

Scalar flavored DM is somewhat more strongly con¬ 
strained than its fermionic analog because of the large 
cross section for producing the colored fermionic medi¬ 
ators of SFDM at the LHC: they must typically be at 
the TeV scale or higher, with the possibility for lower 
masses > 500 GeV only if the dark matter is heavy, 
TO 0 ~ 300 GeV. As a result, many flavor-changing pro¬ 
cesses are suppressed even without any hierarchical or 
MFV structure in the new flavor-violating Yukawa cou¬ 
plings. Moreover annihilation processes relying upon me¬ 
diator exchange can only be dominant for heavy DM, 
> 300 GeV. 

Also distinct from fermion FDM, scalar FDM can have 
important interactions through the Higgs portal. These 
tend to dominate in DM annihilation processes, and will 
be generated at one loop by the mediator couplings A^ 
even if absent at tree level. We showed as a result that it 
is unnatural to expect mediator-dominated annihilations 
outside of the heavy DM mass range mentioned above. 
Another novel consequence is that there can be destruc¬ 
tive interference between Higgs and mediator exchange 
for DM scattering on nucleons, allowing for relaxation of 
direct detection constraints relative to models with only 
one kind of interaction. This makes it possible for SFDM 
to be relevant for indirect detection by gamma rays from 
the galactic center or dwarf spheroidals, unlike for mini¬ 
mal scalar DM coupled through the Higgs. 

Low-energy data from AF = 2 meson mixing provide 
some of the most stringent constraints on the dark matter 
couplings, summarized in table [l] The couplings A^- for 
i,j = 1,2 must either be very small, or very close to 
being diagonal. This is in contrast to the large values A^a 
required if DM annihilation into top quarks is significant 
for determining the relic density. 

There is one intriguing exception: B^-Bg mixing. The 
measured values of the magnitude and phase of the mix¬ 


ing amplitude M 12 are consistent with the SM predic¬ 
tions. However, because of large experimental errors or 
theoretical uncertainties, there is ample room for a new- 
physics contribution to M 12 . We find that, for reason¬ 
able values of its couplings, scalar flavored DM can con¬ 
tribute significantly to both AMg (= 2 |Mi 2 |) and the 
GP-violating phase Pg (= ±arg(Mi 2 )). 

Another interesting possibility in the case of light dark 
matter is the apparently flavor-violating top quark decay 
t —>■ c4>t4>c^ which does not rely upon any off-diagonal 
couplings since flavor conservation is invisibly accom¬ 
plished by the dark matter flavors. 

An important caveat to our analysis which could de¬ 
serve further study is the simultaneous presence of medi¬ 
ators that couple to both left- and right-handed quarks, 
and which mix with each other. By excluding this more 
elaborate class of models, we found that all amplitudes 
involving mediator exchange were suppressed by 1 /m^ 
(typically times a quark mass) rather than \lm^. But 
in more complicated models with mediators coupling to 
both chiralities, one could expect much larger amplitudes 
involving mediators, both for DM annihilation and for 
FGNG processes. 


Appendix A: Simulation of y —^ events 

Here we give details of our simulation of the production 
and decay of u-like mediators that decay to DM plus top 
quarks. We generate events in MadGraph, interfaced 
with Pythia-6.4 for showering and hadronization. 
Events are generated with up to one additional jet at 
the matrix-element level and matched using a shower-kT’ 
scheme. Events are then passed to PGS [35] to simulate 
detector response and event reconstruction. The anti- 
k-r algorithm is used for jet reconstruction, with jet size 
parameter R = 0.4. Event selection is performed with 
cuts corresponding to the signal regions defined in ref. 
m- The same A-factor of 1.5 is applied to the cross 
section. We choose discrete values of the couplings so as 
to vary the branching fraction for decays to tops. As the 
branching fraction depends non-trivially on the mediator 
mass in the very low-mass region, we give the branching 
fraction as a function of the mediator mass in fig. for 
the chosen values of the couplings. We allow for coupling 
to all three generations, but the distinction between first- 
and second-generation quarks does not affect the result, 
for a fixed value of the top branching fraction. 
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Figure 14: Branching fraction of mediator decays to top 
quarks, Xu ^4’^ for chosen values of quark couplings. Xq 
denotes couplings to light quarks, Xq = Xu = Ac. 
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